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1 INTRODUCTION 



ABSTRACT 

We consider the simultaneous chemical, photometric, and gaseous thermal energy evo- 
lution of elliptical galaxies. The evolution of chemical abundances in the intracluster 
medium (ICM) is set by the differing timescales for gas ejection, via supernovae (SNe)- 
driven winds, from dwarf, normal, and giant ellipticals, and is monitored concurrently. 
Emphasis is placed upon the influence of, and sensitivity to, the underlying stellar initial 
mass function (IMF), star formation efficiency, supernovae Type la rates, supernovae 
remnant (SNR) dynamics, and the most recent advances in stellar nucleosynthesis. Un- 
like many previous studies, we adhere to a wide range of optical (e.g. colour-metallicity- 
luminosity relationship) and x-ray (e.g. recent ASCA ICM abundance measurements) 
observational constraints. IMFs biased toward high mass stars, at least during the early 
phases of star formation, are implicated in order to satisfy all the observational con- 
straints. 

Key words: galaxies: abundances - galaxies: elliptical - galaxies: evolution - galaxies: 
intergalactic medium 



What began with Larson (1974b) has led to a pr ofusion 
of fo ll ow-up SNe-d r iven galactic wind models 



establis h. The key to understanding this observed correla- 
tion was provided by Mathews & Baker (1971) but not fully 
appreciated until Larson (1974b). 

Mathews & Baker (1971) postulated that much of the 
gas in the interstellar medium (ISM) of ellipticals had been 
strongly heated by supernovae (SNe) and driven out by a hot 
galactic wind once the gas thermal energy exceeded that of 
its gravitational binding energy (at some time tGw), thereby 
bringing to a halt the bulk of active star formation. Larson 
(1974b) noted that the binding energy per unit mass of gas 
is higher in the more massive galaxies, and thus these sys- 
tems would retain their gas for a longer period of time before 
reaching £gw, and thus attain a higher metallicity, consis- 
tent with the observed mass-metallicity relationship. The 
subsequent evolution would then be regulated only by gas 
lost from dying stars. An excellent recent review of galac- 
tic winds, from an observational slant, is given by Bland- 
Hawthorn (1995). 
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ern advances in stellar nucleosynthesis, supernova remnant 
(SNR) shock dynamics and thermal evolution, the role of 
dark matter, and even basic observational constraints (both 
in the stellar component of ellipticals and the intracluster 
medium (ICM) of galaxy clusters) makes a re-examination 
of the classic wind model a timely one. We have developed a 
software package, entitled MEGaW=Metallicity Evolution 
with Galactic Winds, to enable us to study the chemi- 
cal evolution of ellipticals within the framework of Larson 
(1974b) classic galactic wind model. 

Concurrent to the development of sophisticated wind 
models was the flourishing field of spectral and photometric 
evolution of galaxies, the early history being traced in Tins- 
ley's (1980) seminal paper. One need only look to the recent 
models of Bruzual & Chariot (1993), Bressan et al. (1994), 
Worthey (1994) and Einsel et al. (1995), for a few state of 
the art examples. The consideration of chemical evolution 
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together with photometric evolution is imperative as stellar 
evolutionary tracks and spectrophotometric calibrations are 
sensitive to the chemical composition, and elliptical galaxies 
are complex systems with a distribution of stellar popula- 
tions. Models which utilise only solar abundance tracks and 
solar a bundance spectra/colours (e.g 
Volmer mge 1987 



Bruzual 



|Guidcrdoni fc Rocca ^ 
Chariot 1993[ ) fail to explain 



the observed correlation between the i ntegrated colours and 



absolute magnitude of ellipticals (e.g. Faber 1977 



Bressan 



et al. 1994) (as well as, obviously, the metallicity-luminosity 
relationship) . 

Before embarking upon the construction of a photo- 
chemical evolution code suitable for both the elliptical galax- 
ies in question and ICM abundances, it is important to be 
aware of some of the primary observational constraints^: 

• Elliptical CML Relations: First and foremost, for 
the underlying ellipticals, the models must honour 
the observed, present-day colour-met allicity- luminosity 
(CML) relationships (e.g. Faber 1977). i.e. ellipticals, in 
general, show increasing metallicity and redder colours, 
with increasing luminosity. 

• Magnesium Overabundance in Stellar Populations: 

[Mg/Fe] in giant ellipticals is h0.2 -> +0.3, with 

a slight trend toward increasing values with increasing 
l uminosity, albeit wi th a fairly large degree of scatter 
( |Worthey et al. 199^ ). 

• Oxygen Overabundance of the ICM: X-ray observations 
of the hot ICM of four clusters of galaxies shows an 
oxy gen overabundanc e relative to iron of [O/Fe]=+0.4± 
0.3 ( Mushotzky 1994| ). 

• ICM Iron-Luminosity Relation: The ICM iron mass 
(in Mq) increases as a function of a cluster's opti- 
cal luminosity tied up in E+SOs (in Ly ), such that 
M^ M « 0.02Lv +so flArnaud et al. 1992; ). 

• Type la Supernovae Rates in Ellipticals: The present- 
day Type la SN rate in giant ellipticals is i?i a w 0.03 — > 
0.08 SNuQ ( (Turatto et al. 1994 - 



The work described herein is the only study to date 
to draw upon all of the above; the previously mentioned 
models neglect, or fail to satisfy, one or more of these con- 
straints. Photometric self-consistency is the usual constraint 
which suffers, primarily due to a past dearth of non-solar 



calib rations/tracks (altho ugh, see Arimoto 
and 



Bressan et al. 1994) 



Yoshii 1987 



Overcoming this major hurdle 
is only now becoming feasible; specifically, the past two 
years has seen the release of the first wave of fully self- 
consistent isochrones, ranging from ultra-low to sup er metal- 
rich, iso c hrones. These isoch rone compilations (Worthey 
1994 ai d |Bertelli et al. 19*94 ) draw heavily upon Kurucz's 
(1993) extensive grid of model stellar atmospheres. 

Section || introduces the galactic wind framework in 
which MEGaW operates. This includes an outline of the 
relevant chemical evolution equations, as well as a descrip- 
tion of the different parametrisa tion s adopted for each of 
the input ingredients. In Section 3.1 we present a template 
mode l aga ins t whic h future comparisons will be made. Sec- 
tions 3.2.1 to 3.2.€ then systematically explore the sensitiv- 



* Hq = 50 km/s/Mpc assumed throughout, 
t 1 SNu= 1 SN/century/10 10 L B „ . 



ity of the model predictions to the different input ingredi- 
ents. Such an analysis is long overdue, as previous studies 
tended to simply select one preferred combination of param- 
eters without illustrating the ramifications of said selection. 
A summary of the primary conclusions can be found in Sec- 
tion u 



2 THE CLASSIC GALACTIC WIND MODEL 
2.1 Supernova progenitors and rates 



We adopt the "single degenerate" model of Whelan & Iben 



1973 for Type la SNe in which the progenitors are C/O 



white dwarfs in close binary systems accreting material from 
a secondary companion. Type II SNe are presumed to origi- 
nate via the core bounce- induced explosion of single massive 
star (i.e. m > 8 Mq) progenitors. Arnett (1996) provides an 
excellent review of the relevant SN physics. 

As in Matteucci & Tornambe (1987), we have chosen to 
use the SN rate formalisms proposed by Greggio & Renzini 
(1983). Calling mi and 7712 the mass of the primary and 
secondary, respectively (i.e. ttib = mi + 7712), and denoting 
the secondary mass fraction = 7712/mB, we can write the 
Type la SN rate J?i a (t) as: 



f(^{t-T m2 )dfi dm (1) 



Ru(t) = A 



where m = max[2m t , 3.0]. For a given t, m t represents the 
mass of stars currently leaving the main sequence. 4>(rn) is 
the IMF, by mass. The distribution function for the mass 
fraction of secondaries is taken to be proportional to /i 2 , and 
the minimum mass fraction contributing to the SN Type 
l a rate at time t is /i m = m ax[m2(t)/mB, 1 — 0.5tobm/?tib] 
jGreggio fc Renzini 1983 ). A represents the mass fraction 



of the IMF which is tied up in binary systems with total 
masses in the range m,Bm(3.0 Mq) — > ?tibm(16.0 Mq). A is 
fixed a posteriori by ensuring reproduction of the current 
rate of Type la events in ellipticals. 

The Type II SN rate Rn{t) is composed of two terms - 
one for all stars with initial masses greater than 16.0 Mq, 
and one for the fraction (i.e. 1 — ^1) of stars in the mass range 
8.0 — * 16.0 Mq which are not part of binary systems: 



/•"IBM 

Ru(t) = {I -A) 

J max[m[ ,8.0] 

^ 0(m) 



cj){m) 



4>(t - r m )dm + 



1)1 



-i>{t- 



i)dm. 



(2) 



Single low and intermediate mass stars (m < 8.0 Mq) 
are presumed to end their lives as white dwarfs, after passing 
through a thermally pulsing-asymptotic giant bran ch (T P- 



AGB) and p lanetary nebulae (PNe) ejection phase (Renzini 
fc Voli 198l[). 



2.2 Stellar lifetimes 

Stars return gas and metals to the ISM via stellar winds, 
SNe events and envelope ejection, depending upon the ini- 
tial mass and metallicity. The bulk of the ejection occurs 
near the end of the star's lifetime r(m, Z), with stellar evo- 
lutionary theory providing the functional form of t. 
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Figure [l] illustrates three popular compilations of r ver- 
sus inital mass m: AY87 and AG90 adopted the singular 
power-law form from Talbot & Arnett (1971), which predicts 
excessively long lifetimes for stars in the ~ 1 — »~ 50 Mq 
range; MT87 used Glisten & Mezger's (1982) parametriza- 
tion of the Alcock & Paczynski (1978) and Becker (1981) 

stellar models^]. The best approach available today is 
to adopt the metallicity-dependent lifetimes provided by 
newer, more extensi ve, grids of stellar evolution tracks (e.g. 
t he Padova Group - [Fagotto et al. 1994 or the Geneva Group 
- Schaller et al. 1992). Unless otherwise noted, we use the 
lifetimes from Schaller et al. »s (1992) Z=0.02 grid, although 
for lower metallicity systems, their Z=0.001 predictions are 
adopted. 

Varying the s electio n of stellar lifetime formalism is con- 
sidered in Section 



3.2.1 



Figure 1. Stellar lifetimes r as a function of initial mass m. 
TA71=Talbot & Arnett (1971); GM82=Giisten & Mezger (1992); 
SSMM92=Schaller et al. (1992). 



2.3 Chemical evolution 

Models for chemical evolution follow abundance changes in 
the ISM of a region, and the resulting abundance distribu- 
tions in stars. A detailed derivation of the fundamental equa- 
tions can be found in Talbot & Arnett (1971) and Tinsley 
(1980). We have chosen to couple these equations (appro- 
priate for the one-zone model with instantaneous mixing of 



t Because Alcock & Paczynski (1978) and Becker (1981) did not 
extend their stellar evolution studies below m = 2 Mq , nor con- 
sider metallicities besides Z=0.03, Glisten & Mezger (1982) were 
forced to extrapolate into regimes the models were not designed 
to replicate. This results in considerably shorter lifetimes for solar 
and sub-solar mass stars, in comparison with present-day models. 



gas and relaxation of the instantaneous recycling approxi- 
mation) with the classic galactic wind model formalism of 
Matteucci & Tornambe (1987). 

The mass of the ISM (i.e. gas mass M g (t)) changes 
through star formation and ejection: 



dM g (t) 
At 



-i>(t) + E{t), 



(3) 



where ip(t) is the mass of gas being converted into stars per 
unit time: 



1>(t) = uM 8 (t), 

and E(t) is the total ejection rate of gas from all stars: 



(4) 



E(t) = 



A 



4>{m)^>(t — T m )R(m)dm+ 




f(n)i/){t - T m2 )i?(m 2 )d/i >dm + 



{I -A) 



(f>(rn)t[>(t — T m )R{m)Am ■ 



<f>{m)ij}(t — T m )R(m)dm. 



(5) 



The star formation rate denoted by equation ^ is equiv- 
alent to assuming a Schmidt (1959) Law with exponent 
one. The respective integration lower limits to equation ^ 
are m 1 = max[mi,mL], m 2 = max[mt,mBm], and m 3 = 
max[mt, ttibm]. R(m) = [1 — w m /m] represents the frac- 
tional mass of a star of initial mass m and remnant mass 
Wm, ejected back into the ISM after its stellar lifetime r. 
The lifetime r and remnant mass w m both depend upon a 
star's metallicity Z. Unless otherwise noted, we use the rem- 
nant mass formalism of Prantzos et al. (1993). Following the 
majority of the earlier studies, we take the star formation 
rate if) during the post-galactic wind phase to be zero^J 

The first integral on the right-hand side of equation ^ 
represents the rate at which gas is being returned to the ISM 
at time t from single low mass (m < 8 Mq) stars ending their 
lives as white dwarfs. Integral two is the gas mass return rate 
from stars in binary systems which end their lives as Type 
la SNe. Integral three is the rate of gas mass ejection at time 
t from single stars in the mass range 3.0 — » 16.0 Mq. For 
m < 8.0 Mq, these single stars end as white dwarfs, whereas 
for 8.0 < m< 16.0 M they end as Type II SNe. The final 
integral is the gas ejection rate from single massive stars (i.e. 
16.0 M — > mo) which end their lives as Type II SNe. 

The equation governing the evolution of the mass of 
metals M z (t) = M g (t)Z(t) in the ISM gas is 



3 For the broad-band colours to which the current code is aimed 
at reproducing (e.g. B-V, V-K, etc), relaxing such a restriction 
does not dramatically alter the model results; Arimoto (1989), 
Ferrini & Poggianti (1993), and Elbaz et al. (1995), each provide 
evidence for the limited importance of post-tQw star formation. 
On the other hand, such an assumption is most likely not an 
ideal one when o ne is more concerned with replicating specific 
line indices (e.g. Wbrthey et al. 1995| ); because we are still in 
the process of extending MEGaW's functionality to include full 
spectral synthesis, we have chosen the more conservative ip(t > 
tow) = route, for the time being. 
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dMz(t) 
dt 



= -z(t)m + Mt), 



(6) 



where Ez(t) is the total ejection rate of new and old metals 
(i.e. processed and unprocessed, respectively). The appropri- 
ate form for E%{t) can be reco vere d by replacing mR(m) in 
equation fA with rn^ (Section 2.4 ) , the total mass of metals 



ejected from a star of initial mass m and initial metallicity 
Z(i - r m ). 

As a final sanity check on our numerical solutions of 
equations |§] and |(| we were fortunate to find several experts 
in the field willing to run well-defined standard models with 
which to compare with our o wn code - speci fically, Francesca 
Matteucci (SISSA/Trieste - |Matteucci 1992) ), Frank Timmes 
(Chicago - Timmes et al. 1995[ ) and Leticia Carigi (CIDA ■ 



Carigi 199^ ). Disregarding minor differences incurred by the 
various approaches to nucleosynthesis in the highly uncer- 
tain m ~ 8 — > 11 M range, and assumptions regarding the 
fate of unprocessed metals in low mass star ejecta, the results 
of the intercomparison were more than satisfactory, and put 
this author's mind at ease. The biggest single discrepancy 
occurred in the magnesium evolution, with Matteucci's code 
predicting lower values, a result which was anticipated, and 
simply reflects a difference in the adopted Mg yield between 
the two codes, and not a numerical problem. We return to 
this difference in subsequent sections. 



2.4 Stellar nucleosynthesis yields 

The key ingredient in any chemical evolution code will obvi- 
ously be the adopted nucleosynthetic stellar yields (i.e. the 
mass of element i ejected by a star of initial mass m dur- 
ing the course of its lifetime). These yields, as provided by 
practitioners of stellar evolution and nucleosynthesis, are in 
general parametrised in terms of a progenitor's initial mass 
m, metallicity Z, and "environment" (i.e. single or binary 
system) . 

Type II SNe are responsible for producing the bulk 
of heavy elements; near the end of the progenitor's (Sec- 
tion 2.1) lifetime, during the carbon-burning phase, a He- 
exhausted core (primarily iron-peak nuclei) becomes effec- 
tively isolated from the rest of the star. The subsequent 



core collapse leads to ~ 10° erg being deposited into the 
overlying mantle (neutrino energy), the final result being a 
compact remnant (neutron star or black hole) and an ejected 
envelope, enriched in metals. 

The last few years have seen an explosion of interest in 
massive star (i.e. m > 10 Mq) evolution, the result being the 
newly (or soon-to-be) available metallicity- dependent yields 
of Maeder (1992, hereafter M92), Woosley & Weaver (1995, 
hereafter WW95), and Langer & Henkel (1995, hereafter 
LH95). Improved solar metallicity tabulations from Arnett 
(1991, hereafter A91) and Thielemann et al. (1995, hereafter 
TNH95) have also come on-line. Each of these five options 
are available in MEGaW. 

Each of the compilations start at the approximate lower 
mass cut-off for Type II SNe (i.e. ~ 10 -» 12 M Q ), but have 
very different upper mass limits. M92 goes as high as 120 
M for both Z=0.001 and Z=0.020; A91 goes to 85 Mq, 
but only for Z=0.020; LH95 have an upper limit of 50 M 
for their models, with Z=0.002 and Z=0.020; WW95 have 
the best metallicity coverage (Z=0.0000, 0.00002, 0.0002, 



0.0020, and 0.0200), with an upper limit of 40 M Q ; TNH95 
only go as high as 25 Mq, and only for Z=0.02. We mention 
in passing that WW95 consider three different models for 
m > 30 M , their so-called A, B, and C models. These 
differ in the amount of energy imparted by the piston in 
their models at explosion initiation. Following Timmes et al. 
(1995), we use the "B" models in this mass regime. 

While each source provides the yields for H, He, Z, C, 
and O, the remaining elements under consideration here are 
covered by some, but not by others. N is listed by all ex- 
cept A91, while Mg, Si, and Fe were not part of the M92 or 
LH95 tables. For the sake of self-consistency with the stellar 
evolution models, we do not attempt to "mix-and-match" in 
order to "fill-in" those elements missing in one compilation 
with those found in another, as the input physics between 
models is usually incompatible from one to another. The 
mix-and-match approach is a necessary one for those using 
the M92 yields, as is seen in Carigi (1994) and Giovagnoli & 
Tosi (1995) - e.g. for the evolution of Fe, the former adopts 
a yield of 0.075 Mq for all Type II SNe, irrespective of pro- 
genitor mass or metallicity, while the latter adopts the Fe 
yields of Arnett (1978). 

A detailed comparison of the five Type II SNe yield 
options will be published elsewhere (Gibson 1997, in prepa- 
ration). This will include a comprehensive outline of the dif- 
ferent assumptions regarding the input physics. We refer the 
reader to each of the relevant references for tabulated values 
of mf m (equation ||), although we do draw attention here 
to several interesting points which can be made from a cur- 
sory glance at some of the predicted yield ratios from said 
compilations. 

• [O/Fe]: Mushotzky (1994) has demonstrated that the 
intracluster media (ICMs) of the four clusters in their 
ASCA analysis, is overabundant in oxygen when com- 
pared to iron, by factors of ~ 2 — > 5, with respect to the 
solar ratio. This observation is reflected by the shaded 
region of Figure |j| Contrast this constraint with the 
predicted [O/Fe] yield ratios from the three compila- 
tions which include iron (i.e. WW95,TNH95,A91) - it 
is readily apparent that galactic ejecta enriched in the 
byproducts of Type II SNe of mass m > 20 Mq are 

the necessary primary^] contributor of the bulk of these 
heavy elements in the ICM. The yield ratios adopted 
in Matteucci's (1992) chemical evolution code are in- 
dicated by the dotted line in Figure ^|JJ] For m > 30 
Mq, A91's different mass cut results in more iron being 
ejected for a given initial mass (Gibson 1997, in prepa- 
ration). TNH95's significantly coarser grid makes a de- 

^ Because the other primary source of h eavy element enrichmen t 
-i.e. Type la SNe -have [0/Fe]« -1.53 ([Thielemann et al. 199^), 



they can only act to drive [O/Fe] downward from the observed 
overabundance. As Type la SNe tend to dominate the en richment 
hist ory beyond t > 1 Gyr (Matteucci & Greggio 1986; Section 



3.1), this ICM observation provides indirect evidence that early, 
Type II SNe-driven winds are the dominant mechanism at play. 
II Matteucci (1992) assumed a constant oxygen yield of ~ 0.48 
Mq for all stars of initial mass 9 < m < 15 Mq, suppressing the 
ejection of all other heavy elements. On the other hand, following 
Iwamoto et al. (1994), we have chosen to suppress the ejection of 
all newly synthesised metals in the 8 — > 11 Mq range. 
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tailed comparison difficult, and simple extrapolation to 
higher masses was assumed. 
• [Mg/Fe]: The stellar populations of elliptical galaxies 
posess a ~ 50 — > 100% overabundance of magnesium-to- 
iron, in comparison with the solar ratio (see the shaded 
region of Figure ^ - Worthey et al. 1992). Consid- 
ering the oxygen overabundance relative to iron seen 
in cluster ICMs, with the fact that the bulk of star 
formation in the contributing ellipticals ceases subse- 
quent to the galactic wind epoch, this should not be 
entirely unexpected]^] Significant differences in the be- 
haviour of [Mg/Fe] appear to exist in the stellar mod- 
els in question. Factors of three to ten between A91, 
WW95, and TNH95, at a given mass, are apparent. The 
WW95 curve lies consistently below A91 and TNH95 for 
ra < 35 M . These obvious differences must be borne 
in mind when attempting to replicate Worthey et al. 's 
(1992) observation of [Mg/Fe]^ +0.2 -> +0.3 in ellip- 
ticals. The dotted line in Figure ^ shows the adopted 
[Mg/Fe] in Matteucci's (1992) code.[T] 

The yields for single low and intermediate mass stars 
(i.e. m < 8 Mq) are taken from Renzini & Voli (1981, here- 
after RV81). Stars in this mass range are seen to be impor- 
tant contributors to He, C, and N. These lower mass stars 
are not capable of igniting C in their cores (their masses are 
too low), and end their lives as He or C/O white dwarfs, 
after passing through several dredge-up phases, a thermally 
pulsing phase on the upper AGB, and finally a PN ejection. 
We generally use RV81's preferred models with the Reimers' 



(1975) tnass loss parameter on the red giant branch (KGB) 
77 = 0.33, and the ratio ot the mixing length to the pressure 
scale height a — 1.5]^| 

Stars with initial masses between ~ 8 and ~ 11 Mq 
undergo non-degenerate C-burning and develop O+Ne+Mg 
cores. Whether these highly enriched cores contribute their 
newly synthesised metals to the ISM via thermonuclear ex- 
plosion, or else undergo core collapse, trapping the yields 



** Ty pe la SNc contribute M s and Fe in the ratio [Mg/Fe] « 
-1.63 (Thielemann et al. 1993). Similar to the ICM [O/Fe], this 



lends support to the notion that a Type II SNe-driven wind, as op- 
posed to a later Type la-driven one, is dominating the evolution, 
tt Matteucci's [Mg/Fe] (Figure ^) is fairly flat, and, more impor- 
tantly, negative for all m < 35 Mq. This would appear to be due 
to a lower adopted magnesium yield in her code, which was based 
upon the state-of-the-art during the mid-1980s. Recalling that the 
[Mg/Fe] of the stellar component of giant ellipticals is typically 
+0.2 — > +0.3, one can anticipate that fairly extreme IMFs, with 
a substantial m > 35 Mq component, would be necessary with 
her code, in order to recover such super-solar Mg overabundances. 
This is indeed what Matteucci (1994) found. 

tt Specifically, Tables 3d and 3i (Z=0.020 and 0.004, respec- 
tively) from RV81 for m > 4 Mq , and Tables 3a and 3h (Z=0.020 
and 0.004, respectively) for m < 4 Mq, were adopted. The RV81 
models are admittedly outdated, and our code has subsequently 
been updated to accommodate RV81, as well as the newer Marigo 
et al. (1996) and van den Hoek & Groenewegen (1997) compila- 
tions. None of the results presented herein depend upon the se- 
lection of low mass stellar yields, and we have simply retained 
RV81 for the discussion which follows. A companion paper (Gib- 
son & Mould 1997) contrasts these compilations, as pertaining to 
chemical evolution of the Galactic halo. 



Figure 2. Comparison of the Type II SN e ej ecta [O/Fe] from 
three of the sources discussed in Section [ZA\ The shaded re- 
gion represents the observed ICM [O/Fe] (Mushotzky 1994). Also 
shown is the yield ratio adopted in Matteucci's (1992) chemical 
evolution code. 



in th e remnant neutron star, is still unclear (e.g. 



Matteucci 



1991). Much depends upon the assumed electron capture 
onto 24 Mg and 20 Ne, but the controversy does seem to 
be converging toward the picture in which core collapse 
is favoured, and the elemental enrichment is basically re- 
s tricted to that dep osited during the pre-SN stellar wind 
(Iwamoto et al. 1994), and this is what has been assumed in 



MEGaW. For comparison, Matteucci (1992) suppresses all 
newly synthesised metals except oxygen, and Timmes et al. 
(1995) simply interpolate linearly between the highest mass 
in RV81 (i.e. 8 M ) and the lowest mass in WW95 (i.e. ~ 11 

Mq). 

The yields for each Type la event are taken from the 
improved Z=0.02 W7 Model presented in Thielemann et al. 
(1993). Of the ~ 1.36 Mq of metals ejected per event, iron 
is by far the largest contributor (~ O.74M ), with lesser 
amounts being ejected in the form of Si, O, C, and Mg (i.e. 
0.15, 0.14, 0.05, and 0.01 M , respectively. These yields are 
assumed to be metallicity-independent - a reasonable as- 
sumption based upon the similarity of Thielemann et al. 's 
(1993) Z=0.02 and Z=0.00 models. 

Unless otherwise noted, the default yields adopted 
throughout this paper are culled from WW95 (m > 10 Mq , 
RV81 (m < 8 M ), and Thielemann et al. (1993) (for Type 
la SNe). How the predicted chemical evolution is affected by 
a choice of yield tables d ifferen t from that of WW95 is dis- 
cussed further in Section 3.2.2, and in particular, by Gibson 
(1997, in preparation). 



2.5 Thermal and binding energy evolution of 
the ISM 
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2.5.2 ISM binding energy 



Figure 3. Comparison of the Type II SNeejecta [Mg/Fe] from 
three of the sources discussed in Section Eljj The shaded region 
represents the observed elliptical galaxy [Mg/Fe] (Worthey et al. 
1992). Also shown is the yield ratio adopted in Matteucci's (1992) 
chemical evolution code. 



2.5.1 ISM thermal energy 

For gas to be expelled from a galaxy, the thermal energy of 
the gas heated by SN explosions and st ellar mass loss must 
overcome the binding energy of the gas (Larson 1974b) - i.e. 
the galactic wind will start at a certain time tew when 



£th(tc 



-Eth Ia (iGw) + -Eth n (£gw) + -Bth w (tew) 
fig(tGw), (7) 



where Eth(t) is the total thermal energy in the gas at the 
time t, given by the sum of the contribution from SNe Types 
la and II (E t h lB _(t) and Eth u (t)i respectively) and ther- 
malised kinetic energy from mass loss in high mass (m > 12 
Mq) stars (Eth w (t), which is given by equation 5 of Gibson 
1994a). 

The SNe thermal energy components of equation (?] can 
be written (ISaito 1979bh: 



£th Ia (t) = / £th SN (t-Oi?ia (*')<}£' 



(8) 



Eth n (t) 



£ th stq (t - t')R u (t')dt' 



where t' is the SN explosion time, i?i a and Ru are the SNe 
la and II rates (equations |l| and respectively), and e t h SN 
is the equation governing the evolution of the thermal en- 
ergy con tent in the interior of a supernova remnant (SNR) 
(Section 2.6). In the post- wind phase (t > tGw), the lower 
limit on the time integrals in equations M is taken to be 4gw- 



In order to determine the onset of the galactic wind, we need 
to compute the binding energy of the gas as a function of 
time, fig(t). As discussed in Matteucci (1992), Q s (t) is in- 
fluenced by the presence of dark matter, and its distribution 
relative to the luminous (i.e. gas + stars) component. We 
consider two different scenarios: (i) one in which the two 
components trace each other, following the prescription of 
Saito (1979b), and (ii) one in which the luminous compo- 
nent is embedded in a massive, diffuse halo of dark mat- 
ter, adopting the two-component, self-consistent models of 
Bertin et al. (1992). 

Dark matter traces luminous component 

Following Arimoto & Yoshii (1987), Matteucci & Tor- 
nambe (1987), and Angeletti & Giannone (1990), we first 
represent elliptical galaxies as spheroidal, homogeneous sys- 
tems, with a characteristic gravitational radius Rg, and to- 
tal mass Mo = Ml + Mo (i.e. the sum of the luminous and 
dark components). Assuming the virial theorem holds, the 
total binding energy of the system £7g can be written as 



fie = 



GM G 
2R g ' 



(9) 



Following Saito (1979b), the binding energy of the 
gaseous component f2 g (t) alone can be written 



fig(i) = ( 



M e (t) 
M G 



M s (t) 
M G 



(10) 



It is this binding energy which is compared with the gaseous 
thermal energy calculated from equation |^ in order to de- 
termine the time of galactic wind onset tGw- 

Diffuse dark halos 

While the previous section's analysis is suitable for 
models in which the dark matter component is distributed 
similarly to the luminous component, it is not suitable for 
the more generally accepted scenario in which the gas and 
stars are embedded in a massive diffuse halo of dark matter 
(e.g. Bertin et al. 1994, and references therein). To model 
such distributions, we follow the prescription laid out in 
Bertin et al. (1992). In this context, the binding energy of 
gas is expressed as 



fig(t) = fi L (t) + n L D(t), 



(11) 



where the gravitational energy of the gas due to the lumi- 
nous matter J7l (dropping the implicit time dependency) 
is 

GM g M L 



2i? L 



(12) 



and the gravitational energy of the gas due to the dark mat- 
ter JIld is 



GMgMp , 

S'LD = = "LD- 

-TIL 

The interaction integral f2 LD is estimated to be 



1 R-l 
2tt R d 



1 + 1.37 



Rl 
Rd 



(13) 



(14) 



Rd and Md is the radial extent and mass of the dark matter 
halo. Md ^ (2 — > 10)A/l and Rd > 5i?L are favoured in 
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these t wo-component models when applied to large samples whic h contains virtually all of this thermal energy (Cox 



of elliplicals (Saglia et al. 1992) 



Mass-energy-radius relations 

One last piece of the puzzle necessary to solve for the 
gaseous binding energy is obvious upon re-inspection of 
equations [| [l^ and [n] - a relationship between the galaxy 
mass, radius, and binding energy. To this end, we use em- 
pirical relationships derived by Saito (1979a) from fitting 
Michie-King models to the surface brightness distribution 
and line-of-sight velocity dispersions for a wide range of 
nearby spheroids (ranging from globular clusters to giant 
ellipticals). This first gives the J7l-Ml relation: 



fi L = 1.64 x l(f 



10 12 



[erg] 



(15) 



where Ml is in solar masses. Combining equations g and |l5| 
we can write the virial radius of our spheroidal galaxies in 
terms of the mass as 



Rl = 26.0 



Ml 
10 12 



[kpc 



(16) 



Both these relations were derived assuming a Hubble con- 
stant Ho = 75 km/s/Mpc. Note that one should replace the 
subscript "L" with "G", in equations [TB] and |l6|, when adopt- 
ing the scenario in which the luminous and dark components 
are distributed similarly. 

We should stress the uncertainty in using these empiri- 
cal J7l-Ml--Rl relations which are based upon present-day 
properties of spheroidals. Observational constraints of this 
ilk, for primeval galaxies, simply do not exist. For example, 
if the proto-galaxy has not fully collapsed at the epoch of 
galactic winds, then the binding energy at this point mi ght 
be a factor of ~ 2 less t han that predicted by equation |l5j 



(Arimoto & Yoshii 198£); conversely, mass lost at taw (ei- 
ther impulsively, or gradually over a long timescale) would 
imply that the binding energy at ~ £gw might be anywhere 
from ~ 1 — > 4 times greater than that predicted by equa- 
tion This latter scenario was explored by Hills (1980), 
Vader (1987), and Angeletti & Giannone (1991). We recog- 
nise that use of Saito's (1979a) present-day relationships is 
not ideal, but to be conservative, restrict ourselves to them 
regardless" I. 



2.6 Supernova remnant interior thermal 
energy evolution 

Another of the primary ingredients to any SNe-driven galac- 
tic wind model is the assumed evolution of the thermal en- 
ergy made available to the ISM by each SN event. It is the 
hot, dilute gas in the interior of these SN remnants (SNRs) 



S= If the Hl~ Ml~Rl relation is really a factor of four greater 
than that predicted by equation [l5| then for a given value of the 
star efficiency u, tGW will occur up to an order of magnitude later 
than normally expected. This has the undesired effect of driving 
the colours/metallicities too red/high; increasing v appropriately 
to compensate drives *gw down to comparable values encoun- 
tered using the canonical expression. For the analysis which fol- 
lows, this difference is not important. 



1972| ) 

For our modeling, we have considered a number of ther- 
mal evolutionary scenarios for SNRs, each of which has its 
basis in either the classic models of Cox (1972) and Cheva- 
lier (1974) (the "A" models, below), or the more sophis- 
ticated treatment of Cioffi et al. (1988) (the "B" models, 
below) The latter models incorporate additional radiative 
cooling processes and explicit metallicity effects. The rele- 
vant equations governing each of the models have already 
been outlined in Gibson (1994b, 1995). A cursory treatment 
can also be found in Matteucci (1997). We refer the reader 
to these papers for the minutae, and for this paper we simply 
provide the following qualitative listing: 

(i) Model Aq: SNR shells continue to expand and cool 
radiatively ad infinitum (after Cox 1972 and Cheva- 
lier 1974). The SNR interior thermal energy evolves as 



where t c is the shell cooling time. 



£th SN OC (t/t c ) 

This is the classic scenario adopted by virtually all pre- 
vious studies. 

(ii) Model Ai: Parallels Ao until the interior pressure in 
the SNR is reduced to that of the pressure of the am- 
bient ISM, thereafter merging and becoming indistin- 
guishable from the surrounding ISM. No further radia- 
tive cooling of the interior is considered. 

(iii) Model Bo: Shells continue to expand and cool radia- 
tively ad infinitum (after Cioffi et al. 1988). The late- 
time behaviour of £th SN varies as oc (i/ic)" 10 . The ad- 
ditional factor of (t/t c )~ 0A in £th SN is due to Cioffi et al. 
's (1988) inclusion of radiative cooling, neglected in the 
classic Cox (1972) and Chavalier (1974) studies. This is 
the modern analogue to Model Ao. 
Model Bi: Shells halt both their expansion and ra- 
diative cooling at the ISM merging time alluded to in 
Model Ai's description. 

Model B2: Shells halt their expansion but continue to 

at 



(iv) 

(v) 
(vi) 
(vii) 



(t/tc)' 



cool radiatively ad infinitum (i.e. £th SN 
later evolutionary stages). 

Model B3: Parallels Model B2 until reaching a cool- 
ing time beyond which radiative cooling is no longer 
efficient. 

Model B 3 : Parallels Model B3, unless expanding shells 
start coming into contact, and overlapping with, neigh- 
bouring shells, in which case the expansion term is 
dropped earlier. 

The above models can be compared visually by referring 
to Figure 1 of Gibson (1994b). Unless stated otherwise, we 
shall use Model B2 for the late-time evolution of £th SN (i)- 



3 ANALYSIS 

3.1 The "template" models 

In order to examine the sensitivity of model predictions to 
the various input ingredients, we first describe a working 
template of models which, in general, satisfy the observa- 
tional constraints outlined in Section oj. The pertinent facts 
are listed in Table [t]. We draw attention to the star formation 
efficiency v in column 2; for the chosen set of input ingredi- 
ents v, was treated as a free parameter and chosen to ensure 
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that the present-day CML relationships (columns 7 to 10 
were recovered, as reflected by the solid curves in Figures 
and[^. This parallels Arimoto & Yoshii's (1987) treatment of 
v. Following the procedure outlined in Gibson (1996a), the 
V-band luminosity-weighted metallicity [< Z >]v (column 

10) is computed for each model elliptical. 

The photometric evolution was coupled to the chemi- 
cal evolution as outlined in Gibson (1996a). In general, the 
metallicity- dependent isochrones of Worthey (1994,1995) 
were adopted, although when those of Bertelli et al. (1994) 
were used, the distinction is made. 

Our working template was generated using a time- and 
metallicity-independent form of the Salpeter (1955) IMF, 
with mi = 0.2 Mq and mu = 65.0 Mq. The metallicity- 
dependent yields of WW95 were used for Type II SNe. The 
thermal evolution of the ISM was governed by Model B2 
(Section 2.6). Following Matteucci (1992), diffuse dark mat- 
ter halos with mass and radial extent ratios relative to the 
luminous component of ten were used. 

A binary parameter A = 0.03 (equation [I]) was chosen 
a posteriori to ensure that the Type la SN rate (column 

11) was con sistent with that obs erved in the local elliptical 
population (Iuratto ct al. 1994). The stellar population of 
the giant ellipticals in our template have [Mg/Fe]» ~ +0.15, 
which is only marginally lowe r than the ~ +0.2 — > +0.3 
observed (Worthey et al. 1992). 

Finally, the mass of gas, oxygen, and iron (columns 4, 
5, and 6) ejected in the galactic wind at taw (column 3) are 
also provided in Table W. 

Figure 1 of Gibson (1996b) illustrates the evolution of 
elemental abundance of the primary metals for the M g (0) = 
10 12 M model of Table |l[ 

Figure ^ shows the relevant Type II and la SNe rates 
for this same 10 12 Mq template model. As expected, beyond 
^gw, the Type II rate quickly drops to zero as star formation 
has ceased (and Type II progenitors have m > 10 Mq), 

whereas the Type la SNe |^| continue right up until the 
present-day (the value shown at t — 12 Gyr reflectin g the 
observed value in local ellipticals - Furatto et al. 1994). The 



first Type la SNe do not appear until t « 35 Myr (lifetime 
of an 8 Mq star - the most mass ive secondary allowed in 
our SNe la formalism - Section 2.1), hence the delay in i?i a , 
relative to Rn, seen in Figure [i. 



3.2 Influence of the ingredients 

The following six subsections provide at least a cursory ex- 
amination of the influence of each of the major input ingre- 



11 



Recall from Section 2.1 



that in the binary Type la SNe sce- 
nario, ratios of secondary-to-primary mass of ~ 1 are favoured, 
and that the explosion "clock" is set by the lifetime of the sec- 
ondary. Imposing the chosen 3 Mq minimum-system-mass means 
that beyond ~3 Gyr (i.e. the lifetime of a 1.5 Mq star), one be- 
comes restricted to mass ratios less than unity, with the ratio 
getting smaller and smaller as time goes by, and thus the distri- 
bution function /(/i) of equation hi favours less probable values. 
This results in an expected "break" and subsequent "downturn" 
in the Type la SNe rate at precisely the point seen in Figure 
^. This downturn is also seen in Figure 1 of Greggio & Renzini 
(1993). 



Figure 4. Time-dependence of the Type la and II SNe rates for 
the M g (0) = 10 12 M Q , SNe energy Model B 2 , template. The 
galactic wind epoch *gw (and hence cessation of star formation) 
is noted. 



dients. As the parameter space can be large, we will restrict 
the analysis to just a few pertinent examples, from which 
the importance of the relevant parameter can be addressed. 
Of particular importance will be the effect upon few, as this 
sets the temporal extent of the star formation phase in these 
galactic models, and thus influences heavily the predicted 
present-day photo-chemical properties of the resultant stel- 
lar populations. We shall see that the results are sensitive 
to some of the ingredients, and insensitive to others. 



3.2.1 Stellar lifetimes 

Using the least and most massive template models of Ta- 
ble p , we now vary the stellar lifetime formalism in order to 
view its influence upon the predicted present-day photomet- 
ric and chemical properties. Four different forms were tried - 
the two Schaller et al. (1992) forms (Z=0.001 and Z=0.020: 
SSMM92i and SSMM92 2 , respectively); Glisten & Mezger 
(1982, hereafter GM82); and Talbot & Arnett (1971, here- 
after TA71). Each were shown in Figure [l]. Re-examination 
of this figure should lead one to the intuitive conclusion that 
the results, for all but perhaps the TA71 "single power-law" , 
should be relatively insensitive to the chosen r form. 

Table [2] shows in a more quantitative sense that this is 
true. For discussion purposes, and the sake of brevity, let 
us restrict ourselves to the M g (0) = 10 12 Mq model. Very 
little difference is encountered when switching between the 



2.2 



GM82 and SSMM92i, 2 lifetimes - recall from Section 
that the former was based upon older, super-solar metal- 
licity (Z=0.03), stellar models, whereas the latter allows 
for a more modern prescription for stellar evolution theory, 
and two different metallicities (Z=0.001 and Z=0.020). The 
lower metallicity SSMM92i form leads to marginally later 
wind epochs (~ 5% later), which in turn leads to marginally 
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Table 1. Template Models: Nucleosynthetic yield sources are Woosley & Weaver (1995) for high mass stars, Renzini & Voli (1981) 
for single low and intermediate mass stars, and Thielemann et al. (1993) for binary-progenitor SNe Type l a. Lu minosities and colours 
derived from Worthey's (1994,1995) isochrones. SNe remnant thermal energy follows Model B2 of Section 2.6, which is derived from 
CiofB et al. (1988). A single power law initial mass function, by mass, of slope x = 1.35, and lower and upper mass limits of 0.2 Mq and 
65.0 Mq, respectively, was used, v is the astration parameter for the star formation rate in equation ^. The galactic wind time 4qw is 
in units of Gyrs. A binary parameter A = 0.03 was adopted. Initial dark-to-luminous mass and radial extents of ten were chosen. All 
post-tow ejecta was retained. 



M g (0) 

1.0e6 
5.0e7 
1.0e9 
5.0el0 
1.0el2 



188.9 
209.7 
123.1 



0.006 
0.007 
0.016 



3.5e5 
1.2e7 
1.7e8 



m 6 

2.2e2 
3.2e4 
3.0e6 



l.Oel 
1.5e3 
2.0e5 



M v 

-8.21 
-12.61 
-15.89 



B-V V-K [< Z >h 



0.68 
0.69 
0.74 



2.08 
2.12 

2.44 



-2.28 
-1.56 
-0.51 



46.0 0.077 3.4e9 l.le8 9.3e6 -20.15 0.85 3.04 +0.13 
17.3 0.440 1.7el0 4.9e8 7.5e7 -23.45 0.92 3.33 +0.44 



0.03 
0.03 
0.04 
0.04 
0.05 



more iron being ejected (~ 8% more), due to the increased 
importance of Type la SNe at t > 0.3 Gyr. The stellar life- 
times of Type II SNe progenitors for the Z=0.001 SSMM92 
models are slightly longer than those for solar metallicity, by 
typically ~ 10%. The main reason for this is that the lower 
metallicity models burn at a reduced luminosity, thereby 
lengthening their main sequence lifetimes. 

One might be tempted to infer that the galactic wind 
epoch should be earlier for the higher metallicity r assump- 
tion because of the shorter stellar lifetimes, and therefore, 
earlier SN explosion. As can be seen in Table |[ this is indeed 
the case. 

Not surprisingly, the extremely different form for the 
TA71 stellar lifetimes, already encountered in Figure [j], leads 
to more substantial differences in the output models. The 
second-to-last entry in Table |2| shows that for v — 17.3 
Gyr -1 , the TA71 r-formalism leads to earlier galactic wind 
times (~ 25% earlier) than that found with the SSMM92 
one, regardless of metallicity. In fact, the wind occurs early 
enough for the chosen model parameters, that the metallic- 
ity evolution does not continue to a late-enough epoch to 
ensure that the final stellar population's colours are consis- 
tent with those seen today - e.g. the final V-K is ~ 0.25 mag 
too blue. 

We might naively have expected the opposite behaviour, 
as the TA71 lifetimes range anywhere from to ~6 times 
longer than those found by SSMM92, and thus the typical 
Type II SN explosion would occur later, leading one to per- 
haps expect a later wind. This is not what is encountered, 
the primary reason being that the substantially longer life- 
times encountered using TA71 means that the bulk of the 
Type II SNe explosions are delayed by many tens of millions 
of years in comparison with the SSMM92 lifetimes. This de- 
lay means that the bulk of the "TA71" explosions occur in 
a gaseous medium whose density has been depleted by the 
ongoing star formation during the "delay" . A reduced den- 
sity leads to a greater SNR cooling time, and thus a greater 
thermal energy contribution per SN event, which in turn 
leads to an earlier wind time. 

The final entry in Table ^ shows the implications of 
forcing the M g (0) = 10 12 Mq model to coincide with that 
seen locally in giant ellipticals. This was done by fixing all 
the parameters, save the star formation efficiency u, which 
was reduced from 17.3 Gyr -1 to 7.7 Gyr -1 . This leads to a 
wind epoch which is three to four times later, but it does 
mean that the mean photo-chemical properties of the stel- 



lar populations have enough time to evolve to that seen lo- 
cally. Because the ejection phase occurs so late, it should 
not be surprising to note that even with the reduced star 
formation efficiency, the total gas mass ejected is down 
~ 35%, and the ejecta's [O/Fe] is also reduced by ~ 0.15 
dex (again, due to the increased contribution from the iron- 
important, longer-lived, Type la SNe), in comparison with 
the SSMM92 2 model. 

The choice of stellar lifetime does not influence the 
resultant photo-chemical predictions by more than a few 
percent (provided the discrepant singular power-law form 
is avoided). Future versions of MEGaW will incorporate 
a fully self-consistent metallicity-dependent lifetimes, a la 
Bazan & Mathews (1990), but for the problems at hand, 
this does not appear to be a pressing need. 



3.2.2 Stellar nucleosynthesis yields 

Let us now turn our attention to one of the primary in- 
put ingredients - the Type II SNe nucle osyn thesis yields. 
The basic data were introduced in Section 2.4, to which the 
reader is referred for specifics. Table ^ shows how the galac- 
tic wind epoch changes as a function of yield selection for the 
five sources. We only show the results for one initial mass - 
M g (0) = 10 12 M . 

The star formation efficiency parameter v — 17.3 Gyr -1 
for all entries in the table without an asterisk in the first col- 
umn. Recall that this was the value required to ensure that 
the 10 12 Mq model, in conjunction with the WW95 yields 
and the other parameters in the template models (Section 
3T|), ended up with V-K« 3.35 and [< Z >] v ~ +0.4 by the 
present-day. This is reflected by the first entry to Table || 

One thing we note immediately from Table ^ is that for 
a given v = 17.3 Gyr -1 , taw is not particularly sensitive 
to yield compilation, except that using LH95 leads to wind 
times which are ~ 10% earlier than the other four. This 
can be understood by recognising that over a wide range 
of Type II SNe progenitor masses, the LH95 yields are typ- 
ically ~ 30 — > 40% lower than the others. This is due to 
the combination of their inclusion of stellar winds, and their 
semi-convection treatment (Ledoux criterion and minimal 
overshooting - Gibson 1997, in preparation). 

The top panel of Figure ^ illustrates how the reduced Z 
yield favoured by LH95 manifests itself in the ISM evolution 



- Z g is consistently ~ 30 — > 40% lower. The cooling time of 
a SNR shell increases as a function of decreasing metallicity 
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Table 2. Influence of stellar lifetime selection upon the galactic wind time 4gw (and hence, gas, O, a nd F e. ejecti on, as well as resultant 



3.1. TA71= Talbot & Arnett 1971 



photom etric properties), r den otes the life time formalism ado pted (all other ingredients as in Section 
GM82= Gliisten & Mezger 1982; SSMM92= Bchaller et al. 1992j - subscripts '1' and '2' refer to Z=0.001 and Z=0.020 lifetimes, respectively. 
Superscript * implies that v has been adjusted from that shown in Table [l] in order to recover the same present-day photo-chemical 
properties shown in the earlier template table. 



M g (0) 

1.0e6 

1.0e6 

1.0e6* 

1.0el2 

1.0el2 

1.0el2 

1.0el2 

1.0el2* 



SSMM92 2 
SSMM92i 
SSMM92i 
SSMM92 2 
SSMM92i 
GM82 
TA71 
TA71 



tew 

0.006 
0.006 
0.006 
0.440 
0.462 
0.454 
0.345 
1.476 



3.5e5 
3.2e5 
3.4e5 
1.7el0 
1.7el0 
1.7el0 
2.1el0 
l.lelO 



'"o 

2.2e2 
1.5e2 
2.3e2 
4.9e8 
4.9e8 
5.0e8 
6.5e8 
2.5e8 



'"Fe 

l.Oel 
5.1e0 
l.Oel 
7.5e7 
8.1e7 
8.2e7 
7.1e7 
4.6e7 



M v B-V V-K [< Z >^ 



-8.21 
-8.25 
-8.22 
-23.45 
-23.47 
-23.47 
-23.44 
-23.47 



0.68 
0.68 
0.68 
0.92 
0.92 
0.93 
0.86 
0.92 



2.08 
2.08 
2.08 
3.33 
3.33 
3.35 
3.11 
3.34 



-2.28 
-2.49 
-2.27 
+0.44 
+0.44 
+0.44 
+0.35 
+0.48 



Table 3. Influence of stellar nucleosynthesis yields selection upon the galactic wind time *gw (and hence, gas, O, an d F e. ejec tion, as 



well as r esult ant photometric proper ties). mS denotes the Type II SNe yields so urce fall other ingredien ts as i n Section 3.1 . A91= 
199in^I95=|Langer fe Henkel lggjl WW95=|vvoosley fe Weaver 1995j TNH95=|Thielemann et al. 1996|; M92=|Maeder 1992 



Arnett 



Superscript 

* implies that v has been adjusted from that shown in Table |l] in order to recover the same present-day photo-chemical properties shown 
in the earlier template table. 



M g (0) 

1.0el2 

1.0el2 

1.0el2* 

1.0el2 

1.0el2* 

1.0el2 

1.0el2* 

1.0el2 

1.0el2* 



WW95 
TNH95 
TNH95 
A91 
A91 
LH95 
LH95 
M92 
M92 



tew 

0.440 
0.463 
0.161 
0.460 
0.160 
0.406 
0.736 
0.454 
0.159 



1.7el0 
1.6el0 
2.5el0 
1.6el0 
2.5el0 
2.0el0 
1.5el0 
1.6el0 
2.6el0 



m 



4.9e8 
5.5e8 
9.7e8 
4.6e8 
8.0e8 
2.3e8 
1.9e8 
2.8e8 
4.4e8 



7.5e7 
7.4e7 
8.9e7 
7.3e7 
8.2e7 
n/a 
n/a 
n/a 
n/a 



M v B-V V-K [< Z >] 



-23.45 
-23.38 
-23.36 
-23.39 
-23.36 
-23.45 
-23.48 
-23.40 
-23.37 



0.92 
0.94 
0.91 
0.94 
0.91 
0.93 
0.94 
0.96 
0.92 



3.33 
3.48 
3.35 
3.47 
3.35 
3.29 
3.33 
3.45 
3.33 



+0.44 
+0.49 
+0.42 
+0.48 
+0.40 
+0.28 
+0.30 
+0.43 
+0.39 



(Cioffi et al. 1988; Gibson 1995), and thus the lower ISM Z g 
encountered with LH95 means that radiative cooling of each 
SN event is delayed relative to the other yield compilations. 
This greater energy per SN is what is responsible for the 
slightly earlier *gw. We note in passing that the earlier £gw 
means that the mass of gas ejected in the galactic wind is 
~ 20% greater for a given v. 

Not surprisingly, the LH95 v — 17.3 Gyr -1 oxygen 
ejected at £gw is ~ 55% lower than that found using the 
WW95, TNH95, or A91 yields. A similar lower oxygen ejecta 
mass is found with the M92 yields. Both LH95 and M92 have 
reduced oxygen yields due to stellar winds, and this is re- 
flected in the oxygen ejected at 4gw in Table []. Note though 
that while both LH95 and M92 models have lower oxygen 
ejected at £gw, only LH95's has a similarly reduced "Z" 
yield. M92's Z yield is boosted by a greatly enhanced car- 
bon contribution, especially for solar metallicity and masses 
m>25 M . 

Again, just restricting ourselves to the v = 17.3 Gyr -1 
numbers for the time being, we see that the TNH95 oxygen 
ejecta appears to be ~ 10% greater than the other models 
run without stellar winds (i.e. A91 and WW95). This may 
or may not be so - an uncertainty is introduced when using 
the TNH95 yields because the most massive model in their 
compilation is only 25 Mq. We have just extrapolated be- 
yond this last point in ord er to estimate the most massive 
Type II SNe yields (Section 2/4 ) , which may not be optimal. 



Note that the TNH95 oxygen yields closely parallel WW95 
up to m = 25 Mq; anything beyond that is at best a rough 
estimate. This problem will be ever-present with the TNH95 
yields due to its limited mass coverage. 

The bottom panel for Figure ^ shows that the ISM 
[O/Fe] converges toward ~ +0.0 — » +0.1 beyond t ~ 0.3 
Gyr, regardless of yield compilation. The TNH95 curve is 
~ 0.1 dex greater than the A91 and WW95 at £gw, and 
considerably higher for t < 0.03 Gyr. Again, this is due 
primarily to the uncertain extrapolation to m > 25 Mq - 
the TNH95 iron yield extrapolates to very small values, but 
more importantly, the oxygen yield grows very large beyond 
m ~ 35 Mq. This latter extrapolation tends to drive C/O 
in the earliest phases of the evolution to very low values (see 
the middle panel of Figure |^) . 

This middle panel of Figure |E| is a nice illustration of 
how the uncertainties in the stellar evolution can manifest it- 
self. Whereas we saw that [O/Fe] only spans ~ +0.1 — > +0.3 
dex for the three yield compilations with iron for t > 0.1 
Gyr, the C/O evolution is far more sensitive to the exact 
compilation used. C/O is seen to span almost a full dex 
during the same time period. Unlike the others, the WW95 
curve has an initial decline in the C/O evolution due to the 
high mass Z=0.0000 stars which have very high C/O. The 
upturn in C/O beyond t » 0.1 Gyr is due primarily to the 
increased importance of intermediate mass stars undergoing 
multiple dredge-ups with stellar winds and PNe ejection. 
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Figure 5. Evolution of the ISM metallicity Z g , and C/O and 
O/Fe ratios as a function of time for v = 17.3 Gyr -1 . This v 
value ensures the "Woosley & Weaver (1995)" model's present- 
day photo-chemical properties are consistent with those observed 



locally, ivolution ceases at tGW- 



The enormous metallicity-sensitivity in the high mass car- 
bon yields of M92 leads to the steeper slope for times t < 0.1 
Gyr, when compared with the other four yield sources. 

Only the LH95 yields required a lower value for v, albeit 
the adjustment was minor, the colours being reddened by 
only a few hundredths of a magnitude in V-K. The other 
three (besides the template WW95) required a factor of ~ 2 
increase in the star formation efficiency, the result of which 
was a galactic wind occurring ~ 0.3 Gyr earlier. This had 
the desired effect of "bluing" the predicted V-K by ~ 0.15 
mag, in better agreement with the mean of the observations 
shown in Figure ^. 

3.2.3 Thermal and binding energy evolution of the ISM 

There are a number of inter-related aspects to the ISM ener- 
getics (thermal+gravitational binding) which have already 
been touched upon in Section 2.5. We will now look at a 



few examples which illustrate the sensitivity of tGw to the 
assumed dark matter distribution (diffuse halo (DH) or dy- 
namically dominant (DM)), as well as a couple of "hidden" 
uncertainties that plague many galactic wind models, but 
which are not widely appreciated. Arguments pertaining to 
the influence of assumed individual SNR energetics are de- 



layed until Section 3.2.4 

Table ^ shows how tew, and the relevant present- 
day photo-chemical properties, vary with assumed initial 
dark-to-luminous mass and radial extent ratios. Recall that 
Md/Ml = Ru/Rl = 10 for our template models. It is read- 



ily apparent that the results are not sensitive to the dark 
matter content, provided it is distributed diffusely, follow- 
ing the prescription of Bertin et al. (1992). This is not sur- 
prising given the numerical example of Section 2.5, which 



showed the dominance of the luminous-luminous interaction 
term, as well as the results of Matteucci (1992). Obviously, 
for models in which the dark matter is more centrally con- 
densed (e.g. D=10; R=3), the potential is deeper (equation 
pl|), and it consequently takes longer for the ISM thermal 
energy to build up to the necessary level to overcome the 
binding energy, but the difference is not extreme, and could 
be minimised with a marginal increase in v. 

The situation is not particularly different if one makes 
the assumption that the dark matter is distributed similarly 
to the luminous matter, our so-called dynamically dominant 
(DD) model. In this c ase the gase ous binding energy fol- 
lows equation (after [gaito 1979b[ ) . Table | shows how the 
model predictions for the template 10 12 Mq model vary for 
initial dark-to-luminous mass ratios D=0, 3, and 10 for the 
DD model. 

It can be seen that £gw is not unduly influenced by 
the presence of a DD dark matter component. This is at 
odds with the conclusion of Matteucci (1992) who found 
very late wind epochs (t G w > 9 Gyr) for M g (0) = 10 12 M , 
with D=10 (her Model C2J. The source of the confusion 
can be traced to equation hoi - when the dark matter was 
distributed similarly to the luminous matter, Matteucci had 
M g (0) incorrectly in place of Mg in the denominator outside 
the brackets. The proper form is as shown in eq uation [ lC 
(and was also laid down in its proper form by Ferrini & 



Poggianti 1993, following their equation 1). By using the 



luminous mass, as opposed to the gravitational mass, this 
leads to an order of magnitude overestimation of the gaseous 
binding energy (through equation |l^) . 

In fact, as inspection of column 3 of Table ^ shows, tGw 
actually decreases marginally as one goes to higher dark 
matter contents, simply because the corresponding increase 

in SNe energy efficiency [jffl] outweighs the accompanying in- 
crease in the gaseous binding energy (through equation llOl) . 
This behaviour, while perhaps interesting on a "mathemati- 
cal" level, may or may not be wholly relevant on a "physical" 
level. Specifically, at some point, the validity of continually 
increasing Rg with Mg , for a given Ml , must be called into 
question, as the predicted surface brightnesses will be at 
odds with the observations. On top of this, as stressed by 
Bertin et al. (1994), and references therein, the dynamics 
of elliptical galaxies are better explained by assuming the 
presence of a massive dark halo. For all of these reasons, we 
will henceforth be restricting ourselves to the "dark halo" 
formalism of Bertin et al. (1992), in order to model the evo- 
lution of the gaseous binding energy. 



II The increase in SN energy deposition efficiency with increas- 
ing total mass, simply comes about because in order to honour 
equation [l6| as Mq increases, for a given Ml, there must also be 
a corresponding increase in Rq . By spreading the same Ml over 
a larger volume of radius Rq, the resulting decrease in the hy- 
drogen number density means a corresponding increase in both 
the SN shell cooling time and the ISM merging time, through 
equations 9 and 14, respectively, of Gibson (1994b), and thus an 
increase in £th SN > through equation 12 of Gibson 1994b. 
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Table 4. Influence of diffuse dark matter halo distribution upon the galactic wind time tQW (and hence, gas, O, and Fe, ejection, as 
well as res ulta nt photometric properties), for the Af g (0) = 10 12 Mq model. D= Mq/Ml; R= Ru/Rl- All other ingredients as described 
in Section [O . 

D R v t G w m c J m§ m^ c Mv B-V V-K [< Z >]v < [Mg/Fe] >v 

10 10 17.3 0.440 1.7el0 4.9e8 7.5e7 -23.45 0.92 3.33 +0.44 +0.13 

10 3 17.3 0.719 8.5e9 2.1e8 4.4e7 -23.49 0.93 3.37 +0.49 +0.08 

3 10 17.3 0.382 2.2el0 6.6e8 9.2e7 -23.44 0.92 3.32 +0.41 +0.13 

n/a 17.3 0.357 2.5el0 7.7e8 1.0e8 -23.43 0.92 3.31 +0.39 +0.13 



Table 5. Influence of dark matter, distributed similarly to the luminous component, upon the galactic wind time *gw (and hence, gas, 
O, and Fe, ejection, as w ell as resultant photometric properties), for the M g (0) = 10 12 Mq model. D= Mrj/M L . All other ingredients as 
described in Section 3.1 



D 


3 
10 



v taw 

17.3 0.569 

17.3 0.453 

17.3 0.421 



1.2el0 
1.6el0 
1.9el0 



"'o 

3.0e8 
4.6e8 
5.3e8 



5.6e7 
7.2e7 
8.0e7 



M v 


B-V 


V-K 


[< Z >] v 


< [Mg/F e ; 


23.47 


0.93 


3.35 


+0.47 


+0.11 


23.45 


0.92 


3.34 


+0.44 


+0.13 


23.45 


0.92 


3.33 


+0.43 


+0.13 



One last "hidden" factor which can alter the calculated 
value of tGw, is the formulation used for estimating the hy- 
drogen number density no - a prime ingredient in estimating 
the amount of energy per SN event which is made available 
to the ISM for driving a galactic wind, via its role in setting 
the cooling time for an individual SNR, no should be based 
upon the ISM density at the time of the SN explosion (i.e. 
no an explicit function of time). In order to simplify the en- 



of star formation, and the resultant photo-chemical prop- 
erties). Table ^ illustrates how sensitive tGw is to the SNR 
thermal energy evolution model. For brevity, we restrict our- 
selvesto an analysis of the M g (0) = 10 12 Mq model. 

Recalling that Model B2 has been adopted for the tem- 
plate models (Section 3.1), we see that its entry in Table 
H yields tGw = 0.44 Gyr, with the appropriate present-day 
photo-chemical properties. Model Ao is virtually indistin- 
ergy calculations some early model s (e.g. |Arimoto fe Yoshii| guishable from B 2 for this example, which is not surprising 
^atteucci fe Tornambe 1987 ) used no(t) = no(0), for when a graphical comparis on of the differ ent e t h SN models is 



1987 



all time t > 0. This point was first alluded to by Angeletti & 
Giannone (1990), and more recently by Gibson (1996b). The 
latter reference quantifies its effect as pertaining to Arimoto 
& Yoshii's (1987) models, and the reader is directed there 
for more details. 



3.2.4 Supernova remnant interior thermal energy 
evolution 

One of the more interesting aspects of our work to date 
has been the re-examination of the role played by supernova 
remnant (SNR) thermal energy in p owering galactic winds. 
As already alluded to in Section 2.6, the older galactic wind 
models 



g. Arimoto fe Yoshii 1987] ; Vlatteucci & Tornambe 



1987; Bressan et al. 1994) adopt the classic Cox (1972) and 
Chevalier (1974) formalism for the energetics (our so-called 
"A" models). This form has since been supplanted by the 
more sophisticated models of Cioffi et al. (1988) (what we 
term the "B" models), and our work is the first to incor- 
porate these improvements. Not only is the evolution of the 
individual SNR affected by this new formalism, the influence 
of overlapping shells has at least been treated to first-order. 
This turns out to be a crucial point, which was recognised 
in Larson's (1974b) seminal paper, but again, not fully ap- 
preciated in many of the subsequent detailed models, save 
those of Dekel & Silk (1986), Babul & Rees (1992), and Nath 
& Chiba (1995). Gibson (1994b, 1995) has already examined 
some of the implications of the new formulations, and as 
such, only a few important new points will be made here. 

The adopted energetics form can play a vital role in 
setting the galactic wind time (and consequently the end 



examined (see Figure 1 of Gibson 1994b). Model Bo (Cioffi 
et al. 's 1988 direct analog to the classic Model Ao, but 
including radiative cooling of the interior and metallicity ef- 
fects) leads to very late galactic wind times (4gw ~ 7 — » 8 
Gyr, for giant ellipticals), which for the model presented 
here would imply star formation rates of ip « 15 M©/yr at 
r edshifts z ~ .4, apparently at odds with the observations 
( |5andage 1988; ). 

Models Ai and Bi lead to significantly earlier wind 
times as the late-time evolution of each individual remnant 
differs from the continual energy- loss models Ao and B2. For 
the same v, this of course results in colours/metallicities 
which are too blue/low in comparison with the template 
Model B2. Model B^ illustrates the results when v is reduced 
to 2.9 Gyr -1 , in order to recover the proper photo-chemical 
properties. The later tGw helps these properties, but at the 
expense, somewhat, of the predicted stellar [Mg/Fe]. More 
importantly, the reduced v leads to an order of magnitude 
more gas being ejected at Igw, despite its later occurrence. 

Model B3 is perhaps the best repre sentation of an in- 
dividual SNR's evolution ( Gibson 1994b ) , and the predicted 
wind time and properties are intermediate to the extreme 
models Ai and Bi, and the template B2. Of course, SNRs 
do not evolve in isolation; they eventually come into con- 
tact with neighbouring shells, overlap, and subsequent SNe 
explosions can occ ur in the subse quent rarefied bubble (or 
superbubble - e.g. Tomisaka 1992| ) . Model B 3 takes into ac- 
count, roughly, the shell overlap effects. Similarly, Model B 3 * 
uses a reduced v = 7.7 Gyr -1 to better reproduce the photo- 
chemical properties. This model, as well as the Models B^ 
and B2 illustrate that treatment of the individual SNR ener- 
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Table 6. Influence of SNR thermal energy scenario upon the galactic wind time *gw (and hence, gas, O, and Fe, ejection, as well as 
resultant photometric properties), for the M g (0) = 10 12 Mq model — "A" models are based upon Cox (1972) and Chevalier (1974); "B" 



models are based upon Cioffi et al. (1988). See text for details. All other ingredients as described in Section 3.1. Superscript * implies 
that v has been adjusted from that shown in Table hj in order to recover the same present-day photo-chemical properties shown in the 
earlier template table (specifically, (V-K) and [< Z >]v)- 



Model 


V 






A 


17.3 


0.453 


1.6el0 


Ax 


17.3 


0.115 


2.0ell 


Bo 


17.3 


7.355 


6.5e8 


Bi 


17.3 


0.151 


1.3ell 


BJ 


2.9 


0.922 


1.6ell 


Br, 


17.3 


0.440 


1.7el0 


B, 


17.3 


0.270 


4.5el0 


B^ 


17.3 


0.206 


7.7el0 


B^ 


7.7 


0.539 


6.4el0 



M v B-V V-K [< Z >] v < [Mg/Fe] >v 



4.6e8 7.2e7 -23.45 0.92 3.34 +0.44 +0.13 

4.2e9 3.4e8 -23.22 0.87 3.03 -0.09 +0.10 

1.5e7 3.5e6 -23.55 0.94 3.45 +0.45 -0.16 

3.4e9 2.8e8 -23.30 0.89 3.13 +0.05 +0.10 

3.9e9 4.1e8 -23.32 0.94 3.31 +0.12 +0.04 

4.9e8 7.5e7 -23.45 0.92 3.33 +0.44 +0.13 

1.4e9 1.5e8 -23.40 0.91 3.28 +0.30 +0.13 

2.3e9 2.1e8 -23.36 0.90 3.22 +0.19 +0.12 

1.9e9 2.1e8 -23.39 0.94 3.35 +0.27 +0.10 



Table 7. The mass fraction / m tied up in stars more mas- 
sive than 12 Mq in the four primary IMFs under considera- 
tion in MEGaW. A mass range of 0.2 to 65.0 Mq was as- 
sumed. S55=Salpeter (1955); AY87=Arimoto & Yoshii (1987); 
KTG93=Kroupa, Tout & Gilmore (1993); S86=Scalo (1986). 



IMF 

AY87 
S55 
KTG93 
S86 



f m (m > 12M©) 

0.323 
0.123 
0.061 
0.035 



getics is an integral component in predicting the mass of gas 
ejected at taw- Ellipticals with similar final photo-chemical 
properties can differ in this prediction by up to an order of 
magnitude. 



3.2.5 Initial mass function 

We can anticipate the influence of IMF selection by simply 
looking at the mass fraction tied up in Type II SNe pro- 
genitors (i.e. m > 12 Mq), as, for example, these are the 
primary source of a-elements in chemical evolution models, 
as well as the major contributor to the ISM energetics. Ta- 
ble ^ lists the mass fraction of stars m > 12 Mq for each of 
the four primary IMFs under consideration, assuming a to- 
tal mass range of 0.2 — > 65.0 Mq. One can see that an order 
of magnitude exists between the extrema, with the natter 
Arimoto & Yoshii (1987, hereafter AY87) IMF having ~ 3, 
~ 5, and ~ 10 times the mass of Salpeter's (1955, hereafter 
S55), Kroupa et al. 's (1993, hereafter KTG93), and Scalo's 

1986, Vipreaftpr S&fij TTVITg rpgpprtivply InrkpH intn Typp 



( 

II SNe progenitors of initial mass m > 12 Mr. 



In Table H we start with the template 10 12 Mq model 
(from Section |3T|) , but vary in turn the IMF selection be- 
tween the four listed in Table [?]. Parallel sets of models were 
run for two different S Ne thermal energy forms - Models 
Arj and B2 of Section 2.6. We note that the binary pa- 
rameter A (column 11 and equation ^) is a function of the 
IMF chosen, ranging from ~ 0.01 for the flat AY87 IMF 
to ~ 0.15 for the steep S86 IMF. These values ensure that 



the predicted present-d ay Type la SNe rat e (column 12) is 
7?ia « 0.05 ± 0.03 SNu ( fTuratto et al. 1994|) . 

A primary conclusion to be gleaned from inspection of 
Table is that the galactic wind time tGw (column 3) is not 
particularly sensitive to the IMF, provided the star forma- 
tion efficiency v is kept constant (i.e. v = 17.3 Gyr -1 ), with 
£gw ranging from 0.35 — > 0.57 Gyr. 

A secondary point of interest can be inferred from the 
S55 and AY87 entries to the Model B 2 section of Table §. 
One might naively expect that the flatter IMF (AY87) would 
always lead to an earlier wind time because of its greater 
proportion of Type II SNe, whereas we found that taw oc- 
curred ~ 25% later for the flatter AY87 IMF. This somewhat 
surprising behaviour can be traced to the metallicity depen- 
dence of the SN shell cooling time and ISM merging time 
(Section 2.6) in Cioffi et al. 's (1988) evolutionary formal- 
ism. Metallicity terms were not considered by Cox (1972) 
and Chevalier (1974), which is why this behaviour has not 
been encountered in previous models (nor in the Model Ao 
S55 and AY87 entries to Table |^), which have all been based 
upon these earlier supernova models. A detailed analysis of 
this unexpected behaviour is forthcoming, although a pre- 
liminary accounting can be found in Gibson (1995). 

While tGw may not be overly IMF-sensitive, because of 
the different proportion of low mass stars in the IMFs (i.e. 
those which can effectively lock-up and remove gas from 
possible subsequent enrichment), columns 4 through 6 show 
graphically that the predicted ejecta's mass and abundance 
can be. For example, the S86 and AY87 IMFs, using SNR 
Model Ao and v — 17.3 Gyr - , lead to almost identical taw, 
but the latter model predicts ~ 9 times the gas mass at taw, 
and ~ 70 and ~ 11 times the mass of oxygen and iron, re- 
spectively. This last point is an interesting one - recall that 



[Q/F e1 1CM lies in the range ~ +0.1 -» +0.7 flMushotzky 



1994 ) - the models just mentioned lead to ejecta with [O /Fe] 

of ~ +0.05 (the AY87 IMF) and 0.74 (the S86 IMF). 

Obviously we must fold in a cluster luminosity function be- 
fore claiming anything, but the fact that one of the IMFs 
leads to ejecta which is almost a full dex outside the observed 
ICM [O /Fe] should be a clue to the anticipated difficulty in 
replicating the observations using the steeper IMFs, a point 
which we addressed in Gibson & Matteucci (1997). 

This iron overabundance relative to a-elements in the 
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Table 8. Influence of IMF selection upon the galactic wind time tGW (and h ence , gas, O, and Fe, 
photometric properties) for two dif ferent SNR the rmal en ergy formalisms (Section for the Afp-(O) 



ejection, as well as resultant 
10 12 M w mo del. "IMF" refers 

to the mass function source: S55=Balpeter 19551; AY87=|Arimoto & Yoshii 19871 KTG93=|Kroupa et al. 19931; S86=fecalo 198^ A is 



the binary parameter (equation lj) necessary to recover present-day Type la SNe rate. All other ingredients as described in Section 3.1 
Superscript * implies that v has been adjusted from that shown in Tabic [ijin order to recover the appropriate present-day photo-chemical 
properties (specifically, (V-K) and [< Z >]v)- 



IMF 



tGW 



m^ o M v B-V V-K 
SNe Thermal Energy Model B 2 



[<Z>] 



A Ru < [Mg/Fe] > v 



S55* 


17 


3 





440 


1.7el0 


4 


9e8 


7.5e7 


-23.45 


0.92 


3.33 


+0 


44 





030 





05 


+0 


13 


AY87 


17 


3 





563 


3.5el0 


2 


3e9 


3.4e8 


-23.08 


0.95 


3.97 


+0 


77 





012 





05 


+0 


31 


AY87* 


88 


3 





049 


9.6el0 


6 


6e9 


6.2e8 


-23.03 


0.83 


3.38 


+0 


44 





012 





03 


+0 


34 


KTG93 


17 


3 





450 


1.2el0 


1 


7e8 


4.9e7 


-23.90 


0.86 


2.96 


+0 


20 





060 





06 


-0 


10 


KTG93* 


3 


8 


4 


180 


7.0e9 


1 


0e8 


5.5e7 


-23.90 


0.90 


3.17 


+0 


33 





060 





08 


-0 


38 


S86 


17 


3 





353 


9.3e9 


1 


le8 


5.4e7 


-23.98 


0.82 


2.67 


-0 


01 





150 





06 


-0 


21 


















SNe Thermal Energy Model Ao 


















S55* 


17 


3 





453 


1.6el0 


4 


6e8 


7.2e7 


-23.45 


0.92 


3.34 


+0 


44 





030 





05 


+0 


13 


AY87 


17 


3 





445 


5.1el0 


3 


4e9 


4.7e8 


-23.05 


0.95 


3.94 


+0 


73 





012 





05 


+0 


31 


AY87* 


65 


2 





048 


1.6ell 


9 


8e9 


8.0e8 


-22.93 


0.85 


3.39 


+0 


32 





012 





04 


+0 


32 


KTG93 


17 


3 





579 


9.1e9 


1 


le8 


4.3e7 


-23.92 


0.86 


2.98 


+0 


25 





060 





06 


-0 


15 


KTG93* 


7 


7 


2 


418 


5.5e9 


7 


5e7 


4.8e7 


-23.95 


0.89 


3.13 


+0 


31 





060 





07 


-0 


37 


S86 


17 


3 





467 


5.7e9 


5 


le7 


4.3e7 


-23.99 


0.82 


2.69 


+0 


10 





150 





06 


-0 


31 



ejecta of the S86 model is of course due to the decreased 
significance of the Type II SNe in the steeper IMFs. Type 
la SNe play a correspondingly bigger role as the increased 
binary parameter A would suggest. Not only does the wind 
ejecta suffer from this Fe enhancement, but the predicted 
luminosity-weighted [Mg/Fe] of the stellar population has 
a similar proble m. The observed rat io in giant ellipticals is 
~ +0.2 +0.3 dWorthey et al. 1992[ ), in line with the AY87 
IMF predictions of Table || The S86 and KTG93 values are 
in the range ~ —0.2 — > —0.4, well outside the observations, 
another argument against steep IMFs in ellipticals. 

Studying the v — 17.3 Gyr -1 entries in Table we 
can see that CML predictions for the S55 IMF match the 
mean of the observations (as this is inherent in the tem- 



plate models of Section 3.1). The flatter IMF, because of 
its enhanced enrichment history, by the same tGW, has red- 
der and more metal-rich stellar populations, at the present 
time. Conversely, the steeper IMFs are significantly bluer 
and metal poor relative to the mean of the observations. For 
example, the KTG93 v = 17.3 Gyr -1 models are ~ 0.4 mag 
too blue (in V-K) and ~ 0.2 dex metal-deficient, whereas 
the AY87 models are ~ 0.6 mag too red and ~ 0.3 dex too 
rich. 

To remedy this, we show a series of models with v varied 
in order to best replicate the present-day photo-chemical 
properties of the ellipticals. These models are represented 
with a * in column 1 of Table |[ For the flatter AY87 IMF, 
this means increasing v by a factor of four to five, leading to 
a much earlier wind (£gw ~ 0.05 Gyr as opposed to ~ 0.5 
Gyr). With the earlier wind comes ~ 3 times the mass of gas 
ejected at tew, and a more extreme [0/Fe]~ +0.2 — > +0.3. 

An even more important result is found when we look 
at the attempt to match the KTG93 and S86 IMF mod- 
els with the observations. Specifically, it was found to be 
impossible to redden (or conversely, enrich) these models 
to the observed mean of V-K=3.33 ([< Z >] v = +0.44) 
(for the luminosities involved here). Even reducing v by fac- 
tors of four to five could only redden the colours by ~ 0.2 



mag, which is still ~ 0.2 mag blueward of the mean. By this 
point the wind epoch has shifted to t > 4 Gyr, which for 
the assumed cosmology would imply active observable star 
formation (tp > 150 Mq/vt) a t redshifts z < 0.9, contrary 
to observations (Bandage 1986). By this time the [Mg/Fe] of 
the stellar population has been decreased from an already 
untenable ~ —0.1, to ~ —0.4, even further removed from the 
observed overabundance of magnesium relative to iron. We 
haven't shown the S86 predictions as they are even worse 
than the KTG93 ones. 

This is an important result, and one which cannot be 
remedied within the closed-box formalism for chemical evo- 
lution currently adopted in MEGaW. For these steep IMFs, 
the metallicity increases slowly because of the lack of high 
mass stars in the IMF (at least, "slower" than for the S55 
and AY87 IMFs). Because of this lack of super metal rich 
stars (compared to that seen with the flatter IMFs), all the 
stars formed with Z=0.0 (i.e. all those formed prior to the 
first SNe explosions - t < 0.03 Gyr) tend to "pull" the 
colours/metallicities too blue/low.* * * This would seem to 
indicate that this is partly an artifact of the star formation 
formalism. An infall model in which the mass of the gaseous 
component increases with time according to some accretion 
timescale (as opposed to the closed-box model in which all 
the gas is present and available for star formation at t — 0) 
may be an appropriate mecha nism to "save" the steeper 
IMFs (e.g. Tantalo et al. 1995), although the difficulty in 



* * * This is a parallel manifestation to that of the "G-dwarf 
problem" in the Milky Way, and one from which all closed-box 
models of elliptical galaxy evolution suffer (Worthey et al. 1996). 
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reconciling the [Mg/Fe] may still be problematic . [* * J Infall 
models will be investigated in a future paper. 

More complex scenarios involving bimodal epochs of 
star formation and IMFs are a separate issue, and considered 
in Elbaz et al. (1995) and Gibson (1996a). 



3. 2. 6 Star formation rate 

It is apparent from Figures ^ and ^ that there is consider- 
able scatter about the mean of the observed CML relation. 
A good part of this scatt er is intrinsic (i.e. beyond observa- 
tional error - Faber 1977). One obvious source of scatter can 



most likely be traced to the inherently simplistic handling 
of star formation in MEGaW. 

To illustrate the effect of varying the star formation ef- 
ficiency, we ran models with efficiency parameters v (equa- 
tion |ij) arbitrarily scaled up or down by a factor two, as 
compared with the template values for v listed in Table |l| 
How this impacts upon the galactic wind time, the mass and 
abundance of the ejected ISM, and the resultant present-day 
photo-chemical properties of the remaining stellar popula- 
tion, can be seen in Table ^j, and graphically in Figures ^ 
and^ 

For all the models, the increased efficiency parameter 
(i.e. 2v models) leads to an increased SNe rate, which for 
the more massive models (i.e. M g (0) > 10 9 Mq) results in 
earlier galactic winds, with a correspondingly greater mass 
of gas and metals ejected (up to ~ 50% more). For the lower 
mass (i.e. dwarf) models, the wind still occurs earlier, but 
there is less mass ejected because there is a ~ 4 Myr delay 
before the most massive Type II SNe progenitors explode - 
the increased SFR during this "delay" phase means less gas 
is available for expulsion at £gw- The opposite behaviour is 
seen for the decreased efficiency parameter (i.e. 0.5^ mod- 
ds). 

Figures n and M show the scatter in the observational 
planes for the factor of two variations in v. The majority of 
the observed data points fall within the 0.5f and 2u curves. 
Obviously there are many uncertainties, but we do agree 
with Arimoto & Yoshii (1987) that much of the scatter in the 
observed photo-chemical correlations may be attributable to 
intrinsic scatter in the star formation efficiency 



4 SUMMARY 

We have HpsrrihpH 



in snmp Hptail trip first vprsinn r»f 



MEGa lW. our coupled photometric, chemical, and ISM 
thermal evolution code, based upon the classic framework 
of Larson (1974b). We have attempted to step through, one 
by one, each of the primary ingredients, in order to demon- 
strate how the timeframe for bulk star formation cessation 
(i.e. £gw)j as well as the resultant predicted present-day 
photo-chemical properties, changes, when adjusting any one 

* * * On the other hand, independent supporting evidence 
against steeper-than-Salpeter (1955) IMFs has already been put 
forth by Matteucci & Gibson (1995), Zepf & Silk (1996) and 
Loewenstein & Mushotzky (1996), based upon the significant a- 
element overabundances, with respect to iron, in the hot X-ray 
emitting intracluster gas. 



Figure 6. Solid curve represents the predicted [Z]-My relation for 
the template models of Table [lj The effects of doubling or halving 
the star formation efficiency parameter v, from the template val- 
ues, are also shown. Observational data from Smith (1985) (open 
circles) and Terlevich et al. (1981) (filled circles). 



of these input parameters amongst what appear to be several 
plausible a priori selections. The sensitivity of the results to 
any one parameter has been hidden in previous studies of 
this sort. 

We do not claim to have performed the strictest of sta- 
tistical studies; the analysis shown is meant to be illustra- 
tive, more than anything. What we can conclude at this 
point though is: 

• The present star formation rate formalism (ip oc M g ) 
preclude s the use of IMFs steeper-than-Salpeter (1955) 
(Section |3. 2. 5[). 



Discriminating between early (tGw i$ 0.1 Gyr) and late 
(tew ~ 0.5 Gyr) galactic winds via photo-chemical con- 
straints alone is not possible. For a given assumption 
regarding the efficiency of SNR energy transfer to the 
ISM, we can usually recover the present-day observa- 
tions by varying the star formation efficiency param- 



eter appropriately (Section 3.2.4; Gibson & Matteucci 



1997) 



For non-extreme distributions, the role played by dark 
matter in setting tGw would appear to be less important 
than at first envisaged by Matteucci (1992) (Section 
3 2 3| ). 

An inverse wind phenomenon has been observed 
whereby tew actually increases with increasingly flat- 
ter IMFs. This comes about because of the metallicity 
dependence in the SNR evolution models of Cioffi et al. 
(1988), which tends to reduce the effectiv e ener gy trans- 
ferred per SN event to the ISM (Section |3.2.5|). 



An improved second version of MEGaW is still under 
development. Many enhancements to the basic code are cur- 
rently underway, including most importantly, a full hydro- 
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Table 9. Sensitivity of *GWi an d the resultant present-day photo-chemical properties of ellipticals, to the star formation efficiency 
parameter v. The template model of Table [l] is piven first. Arbitrary scaling of v up and down by factors of two arc listed subsequently. 
All other input ingredients as discussed in Section p.lj . 
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treatment ( Fcrrini & Poggianti 1993 ), the inclusion of a gas 
infall term (Tantalo et al. 1995), and replacing the photo- 
metric evolution module with a full spectral synthesis one. 
A number of other minor enhancements were laid out in 
Gibson (1995). 

Subsequent papers in this series will examine the evo- 
lution of galaxy cluster ICM abundances, damped Lyman-Q 
systems, and stellar yield compilation differences. 
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